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A completely analytical and unified approach to the theory of chirped-pulse oscillators is presented. 
The approach developed is based on the approximate integration of the generalized nonlinear com- 
plex Ginzburg-Landau equation and demonstrates that a chirped-pulse oscillator is controlled by 
only two parameters. It makes it easy to trace spread of the real-world characteristics of both 
solid-state and fiber oscillators operating in the positive dispersion regime. 



1. INTRODUCTION 

High-energy laser oscillators nowadays allow high- 
intensity experiments such as direct gas ionization [l|, 
where the level of intensity must be of the order of 10^'* 
W/cm^. One can expect soon pump-probe diffraction ex- 
periments with electrons and direct high-harmonic gen- 
eration in gases and production of nm-scale structures 
at the surface of transparent materials. Each case calls 
for an intensity of the order of that demonstrated above 
or higher, which means generating tens up to hundreds 
of /LtJ pulses at the fundamental MHz repetition rate of 
an oscillator From the examples above it can be 
seen that high repetition rates are preferable to kHz rates 
(now commercially available) because the signal rates in, 
for example, electron experiments are usually low and 
an improvement factor of 10'^ - 10* due to the higher 
repetition rate of the pulses significantly enhances the 
signal-to-noise ratio. In addition to this physical factor, 
existing kHz systems are more expensive, complex and 
less stable. 

There are a few ways of increasing the oscillator pulse 
energy E, which is a product of the average power and the 
repetition period: by increasing the cavity length and/or 
increasing the power [1, 0, [^, [g, 0, Q ■ The catch is that 
a long-cavity oscillator suffers from instabilities owing to 
nonlinear effects caused by the high pulse peak power Pq. 
The leverage is to stretch a pulse and thereby decrease 
its peak power below the instability threshold. Recent 
progress demonstrating the feasibility of this approach 
has been achieved for Ti:sapphire oscillators operating in 
both the negative- (NDR) 9, 10] and positive-dispersion 
regimes (PDR) [ll|, , for near- infrared Yb-doped solid- 
state oscillators operating in both the NDR @, [l3| and 
the PDR [3, [13]) and for fiber oscillators operating in the 
all-normal disp ersion (ANDi) (i.e. positive dispersion) 
regime Il6j|. 

The fundamental difference between the NDR and 
PDR is that, in the former, the Schrodinger soliton de- 
velops [13] ■ The soliton width T and energy E can be 
expressed as [l8| 

T=Vm/lPo, E^2\p\/-fT, (1) 



where 7 is the self-phase modulation (SPM) coefficient of 
a nonlinear medium (active crystal, fiber, air, etc.), and 
f3 is the net-group-delay-dispersion (GDD) coefficient of 
an oscillator cavity. Since the peak power Pq has to be 
kept lower than the threshold value Pth in order to avoid 
soliton destabilization, one can estimate the maximum 
attainable energy as £' = 2PthT. Energy scaling thus re- 
quires pulse stretching. However, the latter results from 
the substantial GDD growth ]quadratically with energy, 
see Eq. ©]: 

1/3] = E^^/APtt,. (2) 

Hence, the pulse width increases linearly with energy 
(correspondingly, the spectrum narrows with E). As a 
result, i) energy scaling requires a huge negative GDD, 
ii) the soliton obtained has a large width, and iii) it is 
not compressible linearly. 

In contrast to the soliton regime, the pulse is stretched 
in the PDR d9] and its peak power is reduced due to 
chirp tp [20I . [21! ] . The chirp compensates the narrowing 
of the spectrum with energy [20|: 

2A w 87ri/^/7S, (3) 

where A is the spectrum half-width. The advantage is 
that such a pulse [chirped solitary pulse (CSP)] is com- 
pressible linearly down to T w 2/ A (the compression 
factor is « ip). 

As was found, the CSP can be described as a solitary 
pulse solution of the cubic nonlinear complex Ginzburg- 
Landau equation (CGLE) |21| or, more generally, the 
cubic-quintic nonlinear CGLE |20l. [2^. [23l. [24] . This equa- 
tion is thegeneralized form of the master mode-locking 
equation [2lj, [2^, , which provides an adequate de- 
scription of mode-locked oscillators (both fiber and solid- 
state ones). Furthermore, the nonlinear CGLE is used 
in quantum optics, modeling of Bose-Einstein condensa- 
tion, condensed-matter physics, study of non-equilibrium 
phenomena, and nonlinear dynamics, quantum mechan- 
ics of self-organizing dissipative systems, and quantum 



2 



field theory [S^l- Therefore, analysis of the CSP solu- 
tions of nonlinear CGLE is of interest not only from the 
practical but also from the theoretical point of view. 

Since the underlying problem is multiparameter and 
not integrable in the general form, there is no a uni- 
form standpoint on CSPs developing in, for example, the 
ANDi fiber oscillator ^2§\ and the chirped-pulse oscillator 
(CPO) In particular, the physical parameters of os- 
cillators vary greatly and it is not clear a priori whether 
the mechanisms governing the PDR are unified. 

In this work, we propose an approximate method of in- 
tegrating the generalized nonlinear CGLE and show that 
the CSP is its two-parametrical solitary pulse solution. 
As a result, the CSP characteristics are easy to trace on a 
two-dimensional diagram ( "master diagram" ) . Compari- 
son of the PDR parameters demonstrates that the CSPs 
formed in the ANDi fiber oscillator and in the CPO i) lie 
within distinct sectors of the unified master diagram, ii) 
belong to mainly distinct branches of the solution, and 
iii) vary with the parameters in different ways. The vari- 
ation of the main CPO characteristics (spectrum shape 
and width as well as pulse stability) with the PDR pa- 
rameters is analyzed along with the numerical and ex- 
perimental results. A comparison of the models based on 
the different versions of the master equation is made. 

2. CSP SOLUTION OF THE GENERALIZED 
NONLINEAR CGLE 



where z is the propagation distance normalized to the os- 
cillator length (for a ring oscillator model; or to double its 
length for a linear oscillator) , t is the local time [the refer- 
ence frame is co-moving with a solitary pulse solution of 
Eq. ([4])] , a^g — l — ^ is the dimensionless saturated net 
gain {g, I and /i are the saturated gain, non-saturable and 
saturable loss coefficients, respectively), a is the squared 
inverse transmission bandwidth of an oscillator, /3 is the 
GDD coefficient, is the inverse saturation power of the 
self-amplitude modulator, 7 is the net SPM coefficient. 
The normalization of the slowly varying field amplitude 
A{z,t) is chosen such that P{z,t) = \A\ is the instanta- 
neous power. The SAM corresponds to a perfectly sat- 
urable absorber. Such a character of SAM corresponds 
to a semiconductor saturable mirror (SESAM). SESAMs 
are extensively used in the sub- and over-fiJ oscillators 
operating in both PDR [12] and NDR Q. Such a tech- 
nique provides stability and self-staring ability of mode- 
locking and, to date, has no alternative for high-energy 
femtosecond lasers. The main condition, which allows 
describing the SESAM response in the form of (H]), is ex- 
cess of the pulse width (few picoseconds for the PDR) 
over the SESAM relaxation time (wlOO fs) [30]. 

Rescaling z' = zfi, t' = p^j /3, and P' — <;P demon- 
strates that Eq. ^ is in fact three-par ametrical. How- 
ever, as will be demonstrated below, the CSP is the two- 
parametrical solution of Eq. ([4]) . In contrast to Ref. [Ill j 
we do not impose restrictions on the pulse phase (j){t): 



The evolution of the visible, near- and mid-infrared 
electromagnetic fields in an oscillator can be described 
on the basis of the slowly-varying field approximation, 
when the spectral width is much smaller than the carrier 
frequency of the field. The field envelope A is affected 
mainly by i) GDD, ii) SPM (non-dissipative factors), as 
well as iii) saturable gain and linear loss, iv) spectral fil- 
tering, and v) self-amplitude modulation (SAM) (dissi- 
pative factors) pTl [25j . When the effects of higher-order 
dispersion [23| and the field variation along a single oscil- 
lator round-trip 0| are negligible, the oscillator dynam- 
ics obeys the nonlinear CGLE: 



OA ^ , .^.d^A ( /i^ 
— = crA + (a + 2/3) — -- + ^ 



Z7 



\A^A, 



(4) 



A (z, t) ^ ^JP (t) exp (i0 (i) - iqz) , 



(5) 



where q is the phase produced by slip of the carrier phase 
in relation to the slowly varying envelope \J P (t) [17| . 

By analogy with Refs. [13, HI]) substitution of Eq. ([5]) 
in Eq. ^ with the subsequent assumptions a <C /3 (GDD 

prevails over spectral dissipation), and tP^/P ^dt^ <C 1 

(that is, the adiabatic approximation /3 <C T^) results in 



7P(t) = 

dn (t) 



dt 



(3n (^)^ 
n {t) dP (t) 

p (t) dt 



t^^P it) 
1 + <^p (t) 



- an{ty 



(6) 



where (t) = d(j) {t)/dt is the instant frequency. 
Eq. ^ results in 



f3 



dn 

dt 



(A2 - n^) [(A2 - n^) (/i + a) - an^ (a^ 



-7A/5)+(77A/3] 



{A^ -n^ +j/^(3) (A2-3f^2) 
I 



(7) 



The regularity condition dVl/dt < 00 requires that the when n{t)'^ = A^/B. This implies the expression for A^. 
expression in square brackets in Eq. ([7]) equals to zero Simultaneously, the condition P{t) >0 requires the spec- 
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trum truncation so that ri(i)^ < 
arrives at the equations 



q/ (i. Thus, one 



^P(0) = 



fib 
a {E^ 



3 
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2{l + a)-b±VT 

A2 - n^) 



3/3 
2a 
3/i 



= — A^ 



-7A/3 
f 6. 



(8) 



Here a = a/fi, b = aj/iS^n, T = (6 - 2)^+4a (2 + a + b). 
The second equation in ([5]) results from regularization of 
Eq. ([7]) under condition of dVL/dt < oo. 

Positivity of T requires that b < 2 — 2a — 4^/— a, 
with — 1 < a < 0. The condition a < provides 
stabihty against continuum growth, and the inequahty 
a > — 1 means that a cannot exceed the SAM depth. 
The positive ["+"sign in the first equation in Eq. ([H])] 
and negative ("—"sign) branches of the solution ([5]) co- 
incide along the curve P (0) = 3 [2 (1 + a) - 6]/4k = 

3 (y2b -b^l 2b<;. 

Eq. ([5]) demonstrates an existence of two control pa- 
rameters defining the CSP: a and b. Physical meaning of 
the former one is contribution of the saturated net-gain 
a (i.e. saturated gain minus unsaturable loss) regard- 



ing the saturable loss (/i is the maximum saturable loss 
coefficient, i.e. the SAM depth). The /i-parameter is a 
few or fraction of percents for the solid-state CPO, but 
can be substantially larger for the ANDi oscillator (see 
Section [4]). The a-parameter is not free, in fact, because 
it depends on the pulse energy and, thereby, on the b- 
parameter (see below). 

Physical meaning of the 6-parameter is closely related 
to the mechanism of CSP formation [l^, [U [2§| : i) phase 
variation of the propagating chirped pulse can be bal- 
anced in the presence of positive GDD, and ii) spreading 
of the chirped pulse can be balanced by spectral filtering. 
In the cubic nonlinear version of Eq. the CSP chirp 
(-0 >1) is 



(9) 



On the one hand, the peak power Pq relates to the SAM 
saturation power, that is Pq oc i;^^. On the other hand, 
i/'^/T^ c>c A^. Since the CSP spreading is balanced by 
spectral filtering, one can roughly assume A^ cx 1/a. As 
a result of these rough estimations and Eq. ([9]), one can 
see that the ratio b = ^a/f3<;fi = const ^ 1 represents a 
balance of factors providing the CSP existence. 

Returning to Eq. ([8]), its integration results in the 
implicit expression for the CSP profile 



J 



t = 



a,AE{^-A^) { t^^' " - 7] A arctanh (§) + 7^ arctanh (^^^ 



(10) 



where = ii/A^ - ^P{t)/l3 from Eq. Equation 
PH)) allows the spectral chirp to be expressed as 



_ d'^4> [uj) 3/3 A2 - , 



2a (S2-w2)(^2_^2)- 



(11) 



The frequenc; 
pressibility [2. 
chirp is strong 
pulse satellites 
the chirp in the 
border a —Q is 
Ref. [20], the 
frequency lo —0 



dependence of ^ defines the CSP com- 
|3lj: parts of the spectrum where the 
ly frequency-dependent, belong to the 
after pulse compression. The flatness of 
vicinity of a; =0 enhances as the stability 
approached. In contrast to the case of 
chirp is always minimum at the central 



The next step is to assume that 4>{t) is a rapidly vary- 
ing function for the CSP [10, . The stationary-phase 
method of [1^ [s^] allows one to express the spectral 
power from the first of Eqs. 



(12) 



67r/32 A^ ^7/c;/3 
7a 'E? — up' 



e(n' 



2 u?\ 



where 0(x) is the Heaviside function and 11 — min { A, S} 
is the least of A and S. Equation (|12p demonstrates 
that the CSP has the spectrum truncated at ±11 (i.e. 

1^1 <n). 

Integration of Eq. (I12p allows the pulse energy to be 
expressed as 



A 



E . 



6/3^ A 
017 



P {t) dt : 



dw 
p(-)^x 



-A 



(S2 - A2 - 7//3<;) arctanh (I) 
AS 



(13) 
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TABLE I: Relations between dimensional and dimensionless 
quantities 



E E' P^^/'ya^ 

E* E*'(f^ha^ 

p P'A 

A ^'Jll/^ 



TABLE IL Dimensionless equations for the CSP parameters 
(a = o-//i, h = a'y/(3(;ii) 

P' (0) = f [2(l + a)-&±V^] 
A'2 = I [2(l + a)-6±yT] 
r = {b-2f +4a{2 + a + b) 

= I ( A'2 - + b)/ - <^'2) (A'2 - 
p' = 67r(A'2 - + 6) e (A'2 - - c^") 

= 6A' [l - - A'2 - b) arctanh (|;) /a'H'] 



It is clear from Eq. (|13p that the truncation parameter 11 
is equal to A, i.e. A < S. In contrast to the cubic-quintic 
CGLE, whose solution is the truncated Lorentz function 
in the spectral domain [20| . the spectrum in our case is 
parabohc-top (A < S): i) convex if aA^ < 3(/z + a), ii) 
flat-top if aA^ = 3(/x + cr), and iii) concave if aA^ > 

In an oscillator, the cr-(a-)parameter is energy- 
dependent owing to saturation of the gain g. The sim- 
plest law of saturation is g = go/{l + E/ Es) [1^, which 
is valid for the case if the active medium is small in com- 
parison with the confocal length of laser beam {go is the 
gain for a small signal and Eg is the saturation energy). 
Such a law is typical for the fiber and solid-state thin-disk 
oscillators. Reverse relation between the confocal and 
active medium lengthes results in g = g^ ^ ^1 + E/Es 

(3l| . These expressions can be expanded in the vicinity 
of cr = 



(14) 



The £'*-parameter is the energy defined as the averaged 
power of a free-running oscillator multiplied by the cav- 
ity period T,av. E* ^ {g (0) - I - n) EJ{1 + n) for the 
first law of gain saturation (see above and Ref. (26j). 

5 = da/dE\j^^j^, [i.e. 5 = -{I + nf j g {Q) Es] is the 

parameter permitting the gain saturation. For the con- 
sidered law of gain saturation, one has p = E*6 = 
-9 (0) (C - 1)A^ and^^ = g{Q)/{l + fi) ^ g{Q)/l is the 
pump-to-threshold ratio. 

The normalizations presented in Table [T] reduce the 
three-parametrical space of Eq. ^ to a two-parametrical 
one (a, h) for the CSP. The resulting dimensionless equa- 
tions are shown in Table [Til Thus, the CSP is easily trace- 
able on a two-dimensional plane ( "master diagram" ) as in 
the case of the nonlinear cubic-quintic CGLE ,20.. .23.. ,33i] . 



3. CHIRPED-PULSE OSCILLATORS: 
COMPARATIVE ANALYSIS 

The master diagram on the plane (fo, E*) is shown in 
Fig. [TJ The black curve is the stability threshold a =0 
(that \s, E = E*). Above this curve, the CSP does not 
exist (hatched region). Along this curve, the dimension- 
less pulse parameters (see Tables HUT)) are 



< 5 < 2, P'(0) = A'^/b - 3 ( i - i ) , S'2 = 3, 




b arctanh 



6a/1- 



E* 



(15) 

The dashed curve in Fig. [1] shows the border between 
the positive (-I-) and negative (— ) branches of Eq. ([8]). 

The CSPs providing a constant value of the saturated 
net-gain parameter cr correspond to the gray curves in 
Fig. [T] (so-called, isogain curves). The isogain a = is 
the stability threshold (solid black curve), and only three 
nonzero isogain curves are shown (the corresponding val- 
ues of a are superscribed in Fig. [T|) . 

Since cr is a function of E/E*, the isogain curves ex- 

Pin the meaning of the -I— and — branches (see also 
, nil). The -I- branch corresponds to the energy- 
scalable CSP. This means that E grows oc E* along the 
isogain faster than the 6-parameter (6 oc 1//3) decreases 
with the dispersion (3. That is, keeping in such an isogain 
with the energy scaling (note that E* E* / 0^) needs 
a comparatively slow GDD increase. As a result, the 
SPM increases faster than the spectrum degrades with 
the GDD. That is the spectrum broadens (see Fig. [2l 
where 6-decrease along a -I- isogain corresponds to the 
E* -growth in Fig. [T]). The spectra are broadest on the 
stability border. 

For the -I- branch, the spectrum narrows with a (Fig. 
[2] b (X a). This can be explained as a result of the E*- 
decrease, which is necessary for keeping in the isogain 
{E* cx a3/2^* ^ ^pjg^ ^ rpj^g^^ ^j^g gp]^ contribution 

decreases and the spectrum narrows. 

The — branch corresponds to the energy-unscalable 
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~C) 1.0 - 




(6 oc due to growth of the GDD contribution, which 
stretches the pulse when the energy remains ahiiost con- 
stant (Fig. [2]). 

When E* changes weakly along the isogain correspond- 
ing to the — branch, the spectrum broadens with a 
{b oc a) (Fig. [5]). The explanation is that the growth 
of spectral filtering enhances the cutoff of red (bluc)- 
shifted spectral components located on the pulse front 
(tail). The spectral shift at the pulse edges is a conse- 
quence of chirp [l^[2l| . The growth of cutoff shortens the 
CSP and, for a fixed energy, Pq increases. Since Pq (x A^, 
the spectrum broadens. 

One can additionally clarify the division into the -I- 
and — branches by considering the nonlinear cubic limit 
of Eq. (jl]). Such a limit describes a low-energy CPO [2ll |. 
In this case, the exact CSP solution is cx sech(t/T)^+''^, 
and its spectral profile can be expressed through beta 
functions. By the method described in the previous sec- 
tion it is readily found that 



FIG. 1: The master diagram. The solid black curve is the CSP 
stability border. The dashed curve is the border between the 
-I- and — branches of CSP for p/ jj, =17. The gray curves are 
the isogains for a = -0.01, -0.1 and -0.25 (superscribed). 



jP (<) = 13 (1 - tanh^ [At^ (1 + b)/3-f]) , 
Sab 



< 



p{uj) = 



6-2' 

6TTf3 

(TT6K 



(16) 



e (A^ - uj^) 




b 



FIG. 2; The spectral widths in the isogain curves of Eq. (|8]) 
(a = (1), -0.01 (2), -0.1 (3) and -0.25 (4)). The upper (lower) 
branches of the curves correspond to the -I- (—) branches of Eq. 
(ISl): the dashed curve is the border between branches of CSP. 



CSP. This means that b scaling (6 cx 1/(3) weakly affects 
E* {E*' cx E*/(3'^) (Fig. [1]). Thus, the energy remains 
almost constant along this isogain when the GDD scales. 
Certainly, energy scaling is possible as well. However, 
that is not isogain process for this branch of the CSP. For 
the — branch the spectrum narrows with the 6-decrease 



Here ^ = /i<;. Equation (fT6|) demonstrates that i) CSP 
exists only for a <0, ii) its spectral width increases with 
b and \a\. Comparison with Fig. [2] suggests that such a 
behavior corresponds to the — branch. 

On the basis of the model developed, three types 
of CPOs will be considered: i) all-normal-dispersion 
(ANDi) fiber oscillator (lik e that in Ref. [28 |). ii) broad- 
band (e.g. Ti:sapphire [11, llior Cr:YAG[33) and iii) 
narrow-band (e.g. Yb-doped jl3l. [T^ Issjl solid-state os- 
cillators. Representation of these on the master diagram 
affords a means of better understanding and controlling 
of CPO (see also Ref. (H). 



3.1. ANDi fiber oscillator 

In Section [21 the CSP solution of the generalized non- 
linear CGLE was obtained in the limits a <^ (3 and 
/? <C ■ The master parameter b = a.'^ j (3<^pL controlling 
the CSP is defined by the relative but not absolute contri- 
butions from the dissipative and non-dissipative factors 
of the CGLE. This allows a unified standpoint on CPOs 
with parameters which vary within a broad range. 

Figure [3] shows the sector of the master diagram cov- 
ering the ANDi fiber oscillator parameters (see Table [1111 
the estimations for the parameters are based on Refs. 
0, ill). Let us start at point A (Fig. [21 and Table 
IIII[) corresponding to a typical set of ANDi-oscillator pa- 
rameters but with the comparatively small SAM depth 
/i. Although the GDD value is large in comparison with 
that in a CPO, the spectral filter bandwith (25 nm) is 
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FIG. 3: Sector of master diagram (Fig. [T]) representing the 
ANDi fiber oscillators. Parameters corresponding to the A - 
G points are given in TableHni The solid curve is the stability 
border; the dotted and dashed curves are the borders between 
the + and — branches of solution (|H} for p/fj. =3.8 and 0.76, 
respectively (i.e., fi —0.1 and 0.5 for p =0.38, respectively). 



small. As a result, the excess of the ratio (3 /a over that 
for a Ti:sapphire CPO is only tenfold (see below and 
Ref. [13] )■ Simultaneously, an excess of the ratio 7//i<^ 
over that for a Ti:sapphire CPO is tenfold as well. As a 
result, the 6-parameter is ~1 and, dynamically, there is 
no substantial distinction in kind between the ANDi fiber 
and the solid-state CPOs. One difference is that a large 
GDD and a comparatively small E* shift the operational 
point into the — branch region (Fig. [3]). 
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FIG. 4: Numerical (gray circles) and analytical (solid black 
curves) spectra of ANDi fiber oscillators with the parameters 
defined in Table HITl p =0.38. 



TABLE III: Parameters of ANDi oscillator corresponding to 
Fig. 13 E* =12 nJ, 7 =0.014 W"^ The spectrum is centered 
at ~1 /im. 
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A of Fig. |4] shows the spectrum of the numerical so- 
lution of Eq. (Ill) (gray circles) and the analytical profile 
(fT2|) (solid curve) corresponding to point A in Fig. [3l 
The analytical profile reproduces the averaged numerical 
one. As can be seen, the numerical solution is strongly 
perturbed. This effect was first reported in Ref. [37| and 
attributed to excitation of the solitonic internal modes. 
Such modes grow with the GDD and are excited when 
the value (3 /a is large. 

A way of suppressing such a perturbation is to in- 
crease the SAM depth /i (transition 1 from A to B in 
Fig. [3]). The //-growth is not the isogain process, as 



the |a|-parameter increases. As a result, the spectrum 
broadens (Fig. [4l B) according to the analytical model 
(transition from the lower branch of curve 2 to that of 4 
in Fig. O. 

An important factor governing the CPO is spectral fil- 
tering, since the pulse lengthening due to GDD has to 
be compensated by its shortening owing to frequency fil- 
tering of the chirped pulse 21]. In the framework of the 
model under consideration, the filter band narrowing can 
be illustrated by transition 2 to point C in Fig. [31 Such 
a transition results in almost isogain growth of the b- 
parameter. Hence, for the "—"branch of the CSP (Fig. 
^ the spectrum broadens (Fig. d C) [H. The further 
filter band narrowing transforms the pulse into the 
branch of the CSP (Fig. [4l F). Further growth of & oc a 
due to filter band narrowing destabilizes the pulse (cr >0 
at point G). 

The next important factor is the inverse power of the 
loss saturation (the SAM parameter <r) . Its growth (tran- 
sition 3 from C to in Fig. [3]) corresponds to the b- 
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decrease. First, this is not an isogain process, i.e. the \a\- 
parameter increases, which broadens the spectrum (Fig. 
[2|). Second, the spectrum narrows with the 5-decrease for 
the — branch of the CSP (Fig. In our example, point 
D is located in the vicinity of the border between the — 
and + branches of the CSP. In accordance with Figure 
[21 this means that the spectrum broadens owing to the 
|a|-parameter growth (Fig. [H D). D of Fig. [5] demon- 
strates that, unlike the analytical profile, the numerical 
one is distinctly concave. This issue will be discussed in 
Sec. m 

The characteristic of the ANDi fiber oscillator is that it 
is possible to vary the positive GDD within a wide range 
[28l | . The GDD growth decreases the 6-parameter and, as 
a result, narrows the spectrum of the CSP relating to the 
— branch (Figs. [2] and 21 E) 28]. Such a conclusion is 
valid for both isogain and non-isogain variation. In the 
latter case, the |a|-parameter decreases with the GDD 
growth for a fixed E* , which enhances the narrowing of 
the spectrum for the "— " branch (Fig. Solitonic in- 
ternal modes again occur (37j and the spectrum becomes 
perturbed (Fig. HI E). 



3.2. Broadband solid-state CPO 




FIG. 5: Sector of master diagram (Fig. [T]) representing the 
broadband solid-state CPOs. Parameters corresponding to 
points A - D are given in Table IIVI The solid curve is the 
stability border. 



As was pointed out, there are two distinctive differ- 
ences between the ANDi fiber oscillator and the solid- 
state CPO: the former has substantially larger GDD 
and SPM. Nevertheless, both oscillators can be described 
from a unified standpoint because their properties are de- 
fined by only two dimensionless parameters, b and E*. 
From this point of view, the main difference between 
them is that the ANDi oscillator belongs mainly to the 
— branch of the CSP, whereas the CPO belongs to the -|- 
branch. It should be noted that this statement need not 
be considered categorically, because the growth of SAM 
and spectral filtering shifts the operational point of an 
ANDi oscillator into the + branch region (see points D 
and F in Fig. [S]). 

Let us consider a broadband (Ti:sapphire) CPO with 
the parameters presented in Table ITVl The mode-locking 
is provided by SESAM with the inverse saturation power 

which corresponds to a saturation energy flucncc of 100 
/^J/cm^, a relaxation time of 0.5 ps and a mode radius of 
100 /im. 

TABLE IV: Parameters of a Ti:sapphire CPO corresponding 



to Fig. [5] Q =2.5 fs^ 7 =4.55 MW" 
spectrum is centered at ~ 0.8 /^m. 



<;=16 MW" 



The 
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f3 (fs^) 


160 


160 


160 
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0.01 


0.02 


0.02 


0.02 


E* (nJ) 


1120 


1120 


2240 
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The operational point A on the stability border (Fig. 
[5]) corresponds to the minimum GDD and the broadest 



0.15 
CO(fs 



FIG. 6: Numerical (gray circles) and analytical (solid black 
curves) spectra of CPO oscillators with the parameters de- 
fined in Table HVl p=0.n. 



(for a given set of parameters) spectrum (^70 nm, see 
Fig. EA). 

An attempt to increase the energy would destabilize 
the CSP. Therefore, at first it is useful to increase the 
modulation depth (transition 1 to from A to _B in Fig. 
[5]). This is not an isogain process, and therefore the spec- 
trum (Fig. 121 B) does not change substantially in spite 
of the 6-decrease (the 6-decrease broadens the + branch 
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spectrum, but the | a j-growth narrows it, see Fig. [2]). 

The increase of the inverse saturation power <^ (e.g. by 
means of mode reduction) would be useful for subsequent 
energy growth. However, as a rule, the SESAM operates 
in the vicinity of the damage threshold, and so the in- 
growth can be problematic. 

The stability reserve obtained (point B in Fig. [5]) al- 
lows energy scaling (transition 2 from S to C in Fig. O. 
Nevertheless, the twofold energy growth destabilizes the 
CSP in our case, and multipulsing occurs. 

The way out is to increase the GDD (transition 3 from 
C to in Fig. O, which shifts the operational point 
inside the stability range. Since the A ^ D transition 
is an isogain process following the 6-decrease (curve with 
a =0 in Fig. (2), the spectrum broadens (Fig. ^D). 
Again, the spectrum becomes concave. 

It should be noted that, in an experiment, the main 
control parameter used for oscillator stabilization at some 
fixed £'*-level is the GDD value. The GDD growth (cor- 
responding to the 6-decrease) is not an isogain process 
for the H- branch of the CSP and is almost an isogain 
one for the — branch. Figure [2] shows that the spectrum 
narrows in the latter case. For the -I- branch, the \a\- 
parameter increases with GDD, which enhances the CSP 
stability against the continuum growth. Numerical simu- 
lations demonstrate that the effect of the GDD growth on 
the CSP spectrum in the case under consideration is its 
narrowing and stretching of the CSP. Such a conclusion 
agrees with the experiment Figure [7] shows the spec- 
tra from the Ti:sapphire CPO corresponding to the pos- 
itive net GDD growing step by step. The spectra narrow 
with the GDD growth and reshape from concave via flat- 
top to parabolic. Equation reproduces the growth of 
convexity with narrowing of the spectrum. The concave 
spectra appear in the numerical simulations, when the 
spectrum becomes sufficiently broad (Figs. 21 D and [6l 
D). In contrast to Ref. [33|, the present analytical theory 
fails to reproduce this phenomenon (see Section 2]) . 



3.3. Narrowband solid-state CPO 

As was found in Ref. [23|, there is a limit to CPO 
energy growth with resonator lengthening, because the 
modulation depth /i has to increase oc Tcav/Tr (here 
Tr is the gain relaxation time). Since broadband solid- 
state active media such as Ti:sapphire, Cr-doped zinc- 
chalcogenides, and CriYAG have comparatively short 
gain relaxation times (a few microseconds), it is prefer- 
able to use media with a long relaxation time, such as 
Yb-doped crystals. Moreover, the cavity length realized 
so far is already in the MHz-range in terms of repetition 
rate, meaning that the only scaleable parameter remains 
the power. Power scaling is realizable in a Yb-doped 
thin-disk oscillators. For example, the Yb:YAG oscilla- 
tor operating in the NDR has exceeded the 13-/iJ energy 
frontier 8] . Such a regime requires a fair amount of neg- 
ative GDD (w-0.2 ps^ in the case of Ref. ^) and the 



2.0 - 




740 760 780 800 820 840 
wavelength (nm) 



FIG. 7: Experimental spectra from the Ti:sapphire CPO cor- 
responding to the growing positive net GDD (from the black 
via gray to the light gray curves). Pump power 5.5 W; output 
power is of 1.25 W in the mode-locked regime. 

pulse obtained is linearly incompressible. 

It is interesting to consider the prospects of such a 
regime within the PDR. The issue is that a Yb-doped 
medium has a comparatively narrow gainband (a wlOOO 
fs^) and it is not clear a priori that the CPO can operate 
at GDD levels close to a. Nevertheless, Yb-doped CPOs 
have been demonstrated experimentally [H, Il4| . The 
CSPs obtained were compressible down to «450 fs^ and 
the positive net GDD varied within the range ^250 - 
2250 fs^. 

Let us consider a Yb:YAG thin-disk CPO mode-locked 
by SESAM and aiming at an intracavity pulse energy 
level «8G /iJ. At such an energy level, an important fac- 
tor is the SPM due to air filling the resonator 7]. Such 
a factor has to be taken into account. Let the cavity 
length be 15 m, and the averaged mode diameter be 
equal to 2.4 mm. The Yb:YAG-disk thickness is 0.4 mm. 
The SESAM saturation energy fluence is 100 /iJ/cm^, 
its relaxation time is 0.6 ps, and the mode diameter on 
SESAM is 1.2 mm. Other parameters are presented in 
Table El 

TABLE V: Parameters of a Yb:YAG CPO: a =900 fs^ 
? =0.53 MW~\ ^Ji =0.005, E* =80 ^iJ. The spectrum is cen- 
tered at «1 /im. 
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Point A in Fig. |8l corresponds to the narrowband CPO 
operating in the vicinity of the stability border. The 
CPO resonator is filled with air. The CSP belongs to the 
-|- branch; the corresponding spectrum is shown in Fig. 
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FIG. 8: Sector of master diagram (Fig. [T]) representing the 
narrowband solid-state CPOs. Parameters corresponding to 
points A - C are given in Table [V] The solid curve is the 
stability border. 



FIG. 10: Experimental spectrum and autocorrelation func- 
tion (inset) from the Ti:sapphire CPO mode-locked by 
SESAM. Repetition rate is of 60 MHz; intracavity energy is 
of 200 nJ. 
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FIG. 9: Numerical (circles) and analytical (curves) spectra 
of CPO oscillators with the parameters defined in Table FVl 
p=0.Q5. 



[9l A. Its width is «3.5 nm, which allows compressing 
linearly down to w300 fs. Although the ratio /3/a is only 
2.8, the analytical profile provides quite precise fitting of 
the numerical spectrum. 

Violation of the assumption of /3 S> a underlying the 
analytical model results in smoothed spectrum edges. 
Such smoothing increases (Fig. [9l B) when the resonator 
becomes airless and only the active medium nonlinearity 
contributes to the oscillator dynamics. Point B in Fig. [5] 
corresponds to the neighborhood of the analytical stabil- 



ity border, i.e. the almost minimum possible GDD value 
providing the broadest spectrum. Since /? < a, the ana- 
lytical model is not valid, although the pulse is chirped 
and remains almost fourfold compressible (down to ss800 

The net-gainband narrowing (growth of a in compari- 
son with (3) can also result from spectral filtering pro- 
duced by comparatively narrowband SESAM. In this 
case, the a-parameter is defined by the SESAM band- 
width, and the approach of the a-parameter to the /?- 
parameter results in smoothing of the spectrum edges in 
a Ti:sapphire oscillator as well (Fig. [T0|) . 

Spectrum control of the narrowband CPO operating 
in the vicinity of the stability border can be provided by 
inserting a plate (e.g., a sapphire plate) introducing ad- 
ditional SPM (point C in Fig. [51 0.1-cm sapphire plate). 
Figure m C demonstrates that the 6-growth (6 c>c 7) pro- 
vides excellent agreement between the numerical and an- 
alytical solutions. 

It can be concluded that the absolute value of GDD 
required for high-energy pulse stabilization is substan- 
tially lower in the PDR than that in the NDR, so that 
one can avoid the resonator helium filling or vacuumiza- 
tion of the resonator. Nevertheless, vacuumization of the 
oscillator can provide the most direct way to substan- 
tial energy growth. Let us consider an example with an 
intracavity energy of wO.S mJ for a configuration cor- 
responding to point C in Fig. [8l The scaling rules 
E*' = E*ja^/^/^P^, b = aj/(3<;fi (Table [J) demon- 
strate that tenfold reduction of SPM (as a result of res- 
onator vacuumization) and SAM (as a result of, for in- 
stance, mode growth and/or reduction of the saturation 
energy fluence) allow the system to be kept at point C. 
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This guarantees that the dynamics is preserved. How- 
ever, the direct way can be unusable due to Po-growth. 
One can then increase /3 and, simultaneously, decrease 
the modulation depth /i as well as C, 7 (e.g. by means 
of mode growth in a gas-filled resonator) in accordance 
with the scaling rules but keeping Pq below Pth- 



4. BRIEF COMPARISON OF THE MODELS 

Section[2]reports the development of a model of a CPO 
mode-locked by an ideally saturable absorber. Such a 
model is applicable to both the ANDi fiber oscillator (2^ 
and sohd-state CPO mode-locked by SESAM. Simultane- 
ously, there are models of a CPO based on the nonlinear 
cubic-quintic CGLE 0, IH, [H ill, HI]. Such models 
take into account the SAM saturation (i.e. the decrease 
of SAM with power overgrowth), which is important in 
high-power Kerr-lens mode-locked oscillators as well as 
fiber oscillators mode-locked by a polarization modula- 
tor. 

As a result of cubic-quintic SAM, a variety of the non- 
linear regimes appear (25| . In particular, the flat-top 
pulse envelope develops. Since the cubic-quintic SAM 
confines the pulse pe ak p ower, there is the limit for the 
pulse energy growth [23|. It should be noted, that the 
SAM saturation parameter (i.e. parameter defining the 
quintic term in the cubic-quintic SAM) is not truly a 
free parameter in a real- world oscillator. It is closely 
related to the inverse saturation power ^ or to the self- 
focusing power in a Kerr-lens mode-locked oscillator (see 
(20|). In the last case, it is not possible to manipulate 
with SAM and SPM independently because both pro- 
cesses result from the nonlinear refraction in an active 
medium. Therefore such a type of SAM is out of use for 
the micro joule oscillators. 

One has to note, that underestimation of the SAM 
saturation parameter in the cubic-quintic CGLE can lead 
to a huge CSP spectral width and peak power (e.g. see 
[ssjl). Constraint on such a power growth can be provided 
by the gain saturation and, as a result of |a|-growth, the 
CPO operates in the — branch of CPO (in terms of this 
article, see above). Switching to the -I- branch of the CPO 
can be provided by i) ODD decrease, ii) energy growth, 
or iii) spectral filter band narrowing. 

Finally, the low-energy (— branch) sector (as well as 
the large positive GDD sector) of CPO can be described 
by the nonlinear cubic CGLE (i.e. by dissipative gener- 
alization of the nonlinear Schrodinger equation). Such 
a model is useful for a fiber oscillator. For the cubic 
SAM, approach to the + branch of CSP due to the GDD 
decrease can cause the collapse-like instability because 
there is no the peak power confinement due to SAM sat- 
uration. 

An important issue is the appearance of concave spec- 
tra like those in Figs. SI D and [SI D. It is known that, 
for the SAM described by the nonlinear cubic-quintic 
CGLE, the CSP has either parabolic- or finger-top spec- 



tra 0i[23|- However, both numerical simulations [2^.[38j 
and experiment [U, [2^ demonstrate more complicated 
sgectral profiles: concave and convex-concave. In Ref. 
[2J| the analytical concave spectra appear only for the 
cubic-quintic SAM, which enhances the collapse-like in- 
stability (i.e. there is no SAM saturation). In Ref. [33l |. 
the stable analytical CSPs with concave spectra are ob- 
tained for nonzero quintic SPM. For the SAM type con- 
sidered in the present work, the analytical concave spec- 
tral profiles are not admissible because the first equation 
of ([5]) demonstrates that A'^ < 3(1 + a). Hence, the 
existence of concave spectra cannot be explained on the 
basis of the present dissipative soliton models (with the 
exception of the model presented in [s^). 

Figures IH [6] show that the concavity grows with spec- 
trum broadening when 1/(2A)^ tends to a. Simultane- 
ously, there is no concave spectrum, when SPM is sup- 
pressed (Fig. [21 where 7 < ^<;). 

Numerical analysis demonstrates that the appearance 
of concave spectrum is not rooted in the character of 
SAM. Surprisingly, however, the concave spectrum pro- 
file exists even in the case of the nonlinear cubic CGLE 
and develops with the spectrum broadening. Figure [TT] 
shows the spectra corresponding to point B in Fig. [31 

but for the different SAM types: ^i<i\A\'^^ (^1 + <;|A|^^ 

(black curve, see Eq. ([4])) and (nonlinear cubic 

model, gray curve). It can be seen that the absence of 
SAM saturation in the second case leads to spectrum 
broadening due to power gr owth (A^ cx F(0), see Eqs. 
(|8I16[) ). Analytical estimation for A from Eq. (|16p gives 
0.0076 fs~^ vs. numerical value 0.0074 fs^^. As a result 
of the spectrum broadening, the concavity increases. The 
P(t)-profile does not have any anomalies. Thus, growth 
of spectral components at the spectrum edges can be 
treated as a fundamental feature of the regime consid- 
ered (see also Ref. (S^), without obhgatory addressing 
of higher-order dispersions [2§| . 



5. CONCLUSIONS 

Analytical theory of the CSP has been developed. The 
theory treats the CSP as a solitary pulse solution of the 
generalized nonlinear CGLE. The main advantage of the 
theory developed is the possibility of representing the 
CPO parametrical space in the form of a two-dimensional 
master diagram. As a result, the CSP characteristics be- 
come easily traceable. It has been demonstrated that 
both ANDi fiber and chirped-pulse solid-state oscillators 
can be described from a uniform standpoint and repre- 
sented on a unified master diagram. The main difference 
between them is that they realize mainly two different 
branches of the CSP solution. Such branches differ in 
the energy and dispersion scaling rules as well as in the 
behavior of the CSP parameters. 

Comparison with the results of numerical simulations 
has shown that the analytical solution provides a good 
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FIG. 11: Numerical spectral profiles corresponding to the 
model Q (black curve) and its nonlinear cubic limit (gray 
curve). Parameters correspond to the point B in Fig. [31 



approximation of the spectrum shape, which is truncated 
and has a flat or parabolic top. The approximation is 
quite precise even in the case where the net-gainband 
is so narrow that the squared inverse bandwidth verges 
towards the GDD. This provides an adequate description 
of both ANDi fiber and thin-disk narrowband solid-state 
chirped-pulse oscillators. Thus, the theory allows the 
CPO characteristics to be optimized and demonstrates 
the feasibility of at least, sub-mJ energy scaling in a thin- 
disk CPO. 
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